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We introduce equations of motion for spin dynamics in a ferromagnetic Bose-Einstein condensate 
with magnetic dipole-dipole interaction, written using a vector expressing the superfluid velocity and 
a complex scalar describing the magnetization. This simple hydrodynamical description extracts 
the dynamics of spin wave and affords a straightforward approach by which to investigate the 
spin dynamics of the condensate. To demonstrate the advantages of the description, we illustrate 
dynamical instability and magnetic fluctuation preference, which are expressed in analytical forms. 
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I. INTRODUCTION 

One of the salient features of a gaseous Bose-Einstein condensate (BEC) is the internal spin degrees of freedom. In 
spinor BECs, namely, in BECs with internal degrees of freedom, spin and gauge degrees of freedom couple in various 
manners, leading to nontriviai properties of spin waves and topological excitations. For example, ferromametic BECs 
have continuous spin-gauge symmetry, thus the circulation of the superfluid velocity is not quantized whereas 

tin-1 polar BECs and spin-2 cyclic BECs can host fractional vortices due to the discrete spin-gauge symmetry 0- 
, In recent experiments, in situ imaging of transverse magnetization has revealed the real-time dynamics of the 
spontaneous symmetry breaking, spin texture formation, and nucleation of spin vortices openmg up a new 

paradigm for studying the static and dynamic properties of spin textures. 

On the other hand, BECs with magnetic dipole-dipole interaction (MDDI) have also attracted much attention both 
experimentally and theoretically in recent years. The long-range and anisotropic nature of the MDDI is predicted to 
yield exotic phenomena, such as new equilibrium shapes, roton-maxon spectra, supersolid states, and two-dimensional 
solitons [ll|. In particular, when the BEC has spin degrees of freedom, the MDDI is predicted to develop spin textures, 
even when the MDDI is much weaker than the contact interaction |12l-[l^. This work is motivated by experiments 
done by the Berkeley group [sl , where small magnetic domains were observed to develop from a helical spin structure 
in a spin-1 ^^Rb BEC. The method presented in this paper simplifies the spin dynamics in a complicated system of 
spinor dipolar BECs, although we have shown in our previous work that mean- field calculations do not reproduce the 
experimental results [15]. 

In this paper, we propose a new type of hydrodynamic description of a ferromagnetic BEC with MDDI. The 
hydrodynamic equation of spinor BECs has been discussed for both ferromagnetic phases p!6l-[T8| and non-magnetized 
phases [l8|, |T9|. In Ref. [l6|, Takahashi et al. consider the strong MDDI limit by using the classical spin model, i.e., 
by neglecting the spin-gauge coupling. On the other hand, Lamacraft takes into account the spin-gauge coupling by 
introducing the so-called Mermin-Ho relation, and considers the weak MDDI [l^- In these papers, the authors use 
a unit vector to describe the local magnetization in the ferromagnetic phase. Here we use a single complex scalar 
variable instead of a unit vector to describe the local magnetization and treat both the spin-gauge coupling and 
MDDI. This simple description allows a straightforward approach to analyze the spin dynamics of the condensate. 
In order to demonstrate the advantages of our description, we analyze the dynamical instability and magnetization 
fluctuation preference of the BEC with MDDI. 

The rest of the paper is organized as follows. In Sec. [Ill hydrodynamic equations described using the spin density 
vector are derived from the Gross-Pitaevskii (CP) equation. We then rewrite the equations by means of stereographic 
projection for some simple cases: quasi two-dimensional (2D) systems under zero external field and under a strong 
magnetic field. For both zero-field and strong-field cases, the wavevector dependence of the dynamical instability is 
obtained straightforwardly in an analytical form in Sec. Illli In Sec IIV| we also illustrate the magnetic fluctuation 
preference for the unstable modes discussed in Sec. IIIII Conclusions are given in Sec. |Vl 
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II. HYDRODYNAMIC DESCRIPTION 

A. Equations of motion of mass and spins 

We consider a spin-F BEC under a uniform magnetic field B applied in the z direction. The GP equation for the 
spinor dipolar system is given by 

ih^^^m{r,t) = {Ho^pm^qm^)^m 

+ E E E {mn\SMs){SMs\m^n^)K^m^^n^ 

»S'=0,even Ms — — Sn,m',n' — — F 

F 

+^dd E ^M(^M)mn^n, (1) 

fi=x,y,z n= — F 

where ^^(r,t) is the condensate wavefunction for the atoms in the magnetic sublevel m and Hq = — fi^V^/(2M) + 
Utrap{r), with M being the atomic mass and /7trap the spin-independent trapping potential. The linear and quadratic 
Zeeman energies per atom are given by p = QfI^bB and q = (^FMB^)^/^hf , respectively, where is the hyperfine 
^-factor, /iB is the Bohr magneton, and E'hf is the hyperfine energy splitting. The second term on the right-hand side 
of Eq. ([T]) comes from the short-range part of the two-body interaction given by 

2F ,2 

V,{ry) = 5{r-r') ^ ^^sVs. (2) 

5=0, even 

where Vs = Y1^Ms=-s \^^s){SMs\ projects a pair of spin-1 atoms into the state with total spin 5, and as is the 
5-wave scattering length for the corresponding spin channel S. The scattering amplitude for odd S vanishes due to 
Bose symmetrization, and {mn\SMs) in Eq. ([1]) is the Clebsch-Gordan coefficient. The last term on the right-hand 
side of Eq. ([1]) corresponds to the MDDI, where Cdd = Mo(^FMB)^/(47r), with jiQ being the magnetic permeability of 
the vacuum. Here, we define the non-local dipole field by 

b^{r)= fd^r'J2Q^^Ar-r')fAr'), (3) 
where Q^v{r) is the dipole kernel, given in Sec. Ill Ci and 

//i = ^ ^^ ^m^F^rnn^n (4) 
mn 

is the spin density, with F^^y^z being the spin-F matrices. Below we omit the summation symbol: Greek indices that 
appear twice are to be summed over x, and and Roman indices are to be summed over — F, . . . , F. 
From the GP equation, we can immediately derive the mass continuity equation: 



' V • (ntot^mass) = 0, (5) 



dt 
where 

ntot = ^l^^m. (6) 
ntott^mass = ^ [^;^(V^^) - (V^^,)^^] (7) 

are the num ber densit y and superfiuid current, respectively. By introducing a normalized spinor defined by 
^m(^7^) = \/^tot(^, ^)Cm(^, the superfiuid velocity Vmass can be written as 

^ [C(va)-(vC)C,n]. (8) 



2Mi 
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In the absence of the Hnear and quadratic Zeeman effects and the MDDI, the continuity equation of the spin density 
can also be derived from the GP equation as 

^+V.(nto<i„)=0, (9) 

where v^^^^ is the spin superfluid velocity defined by 

<pi„ = ^(J^p)™n[C(VC„) - (VC)Cn]- (10) 

The short-range interaction does not contribute to the equation of motion of spin, since it conserves the total spin of 
two colliding atoms. The detailed calculation is given in Appendix [A] In the presence of the external magnetic field 
along the z direction, the linear Zeeman effect induces a torque term {p/h){z x f)^ on the right-hand side of Eq. ([9]), 
which causes the precession of spins. In a similar manner, the dipole field also induces a torque term (cdd/^)(& x 
On the other hand, the quadratic Zeeman term does not conserve the transverse magnetization and its effect is written 
as {2q/h)e^zv'i^totJ^zv, where Cijk is the Levi-Civita symbol and 

= \CMF^ + F.F^UnCn (11) 

is a nematic tensor. The derivations of these three terms are given in Appendix [A] As a result, we obtain the 
equation of motion of spins in the presence of the MDDI and linear and quadratic Zeeman effects under the external 
field parallel to the z axis: 

^ + V • (ntot<i„) = ^{b X /)^ + |(i X f), + |ntote^..Ar,.. (12) 
Equations (|5|) and ([12]) hold in all phases, independent of scattering length. 



B. Ferromagnetic BEC 

In the following, we consider a ferromagnetic BEC. We assume that the BEC is fully magnetized, |/| = Fritot-, and 
only the direction of the spin density can vary in space. This assumption is valid when the ferromagnetic interaction 
energy is large enough compared with the other spinor interaction energies, MDDI energy, quadratic Zeeman energy, 
and the kinetic energy arising from the spacial variation of the direction of /. The linear Zeeman effect is not 
necessarily weaker than the ferromagnetic interaction, since it merely induces the Larmor precession. For example, 
the short-range interaction ([2]) for a spin-1 BEC can be written as [20] 

{mn\Vs{r,r')\m'n) = 5{r - r') [cQdmn^rn'n' + Ci{F^)mn{F^)m'n'] , (13) 

where cq = 47rfi^(2a2 + ao)/(3M) and Ci = 47rfi^(a2 — ao)/(3M). The ground state is ferromagnetic for ci < 0. The 
above assumption is valid when q <C |ci|ntot7 c-dd ^ |ci|, and the length scale of the spatial spin structure is larger 
than the spin healing length ^sp = ^/ ^2M|ci|ntot- Moreover in the incompressible limit, namely when the spin 
independent interaction (contot for the case of a spin-1 BEC) is much stronger than the ferromagnetic interaction and 
MDDI, the number density ntot is determined regardless of the spin structure and assumed to be stationary. This is 
the case for the spin-1 ^^Rb BEC. 

We then rewrite the equations of motion ([5]) and ([T2]) in terms of a unit vector / = //(Fntot) that describes 
the direction of the spin density and the superfiuid velocity Vmass defined in Eq. ([8]). The order parameter for the 
spin-polarized state in the z direction is given by ("m^ = ^mF- The general order parameter is obtained by performing 
the gauge transformation and Euler rotation as 

^ = g#g-.F.ag-iF,,3g-in7^(0) 

= e»(0-F7)e-iF,ag-in/J^(O)^ (14) 

where a, /? and 7 are Euler angles shown in Fig. [T] and (j) is the overall phase. Due to the spin-gauge symmetry of the 
ferromagnetic BEC, i.e., the equivalence between the phase change and spin rotation 7, distinct configurations of C, 




FIG. 1: Euler rotation of the unit vector /. 



are characterized with a set of parameters (a, f3^(j)' = (j) — Fj). Here, a and /3 denote the direction of / as shown in 
Fig. [TJ Actually, / for the order parameter ([14]) is calculated as 

= C^^^* ((^Fy(5^iF,cxp^-iF,cx^-iFy(5\ M) 
Sm V / rnn 



^sin f3 cos 
F I sin /3 sin a I , 

cos (3 



(15) 



where lZz{(y) and 7ly{P) are the 3x3 matrices describing the rotation about the z axis by a and about the y axis by 
/3, respectively. In a similar manner, we obtain the nematic tensor J\ for the order parameter (p!4|) as 



_F F{2F - 1) . . 



(16) 



where T denotes the transpose and //^^ = ^(m^^* {F/^Fi, + Fi,F^)mnCn^ is the nematic tensor for C,^^\ which is given 
by 



,(0) 



^F/2 
F/2 


^il\t\,Sl^ lilts. 



,0 1 



,0 1. 



Substituting Eq. ([14]) and 



into Eq. ([8]), the superfluid velocity can be written as 



(17) 



(18) 



Vmass = - F{Va) COS/3], 



which satisfies the Mermin-Ho relation l2l|: 



hF 



V X i;„iass = -^e^^xff^C^fu X V/a). 
As we mentioned before, ntot is stationary in the incompressible limit. Thus, Eq. (jSJl becomes 

V • (ntot f mass) = 0. 



(19) 



(20) 



(21) 
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Equations (|2Q|) and ([2T]) are equations for the superfiuid velocity. Next, we consider the equation for the spin superfiuid 
velocity. Making use of Eq. (fT6|) , Eq. ([TOj) can be rewritten in terms of / and i^mass as 



spin 



{V4>')Ff^ - {Va)Af^, - {VpWf^y cosa + 



fix ' 



-Ffv 



hF 



- COS j3 cos a 

- cos P sin 

sin /3 





1— sin a\ 






cos a 


1 




I y 





mass c^y^^fivxfv^ fx- 



(22) 



Substituting Eqs. (fT6|) . ([2T]) . and ([22]) into Eq. ([T2j), we obtain the hydrodynamic equation in terms of / and Vmass as 



dt 



V)/ = -/ X Beff, 



(23) 



with 



eff 



-^(a-V)/-— W 



where a = (Vntot)/^tot and is the unit vector along the external field (e^ = i in this paper). Here we note that 
Eq. (|23|) has the same form as the extended Landau-Lifshitz- Gilbert equation (without damping) which includes the 
adiabatic spin torque term (22| . 



C. Quasi-2D system 



We next consider a quasi-2D system, that is, we consider a BEG confined in a quasi-2D trap whose Thomas-Fermi 
radius in the normal direction to the 2D plane is smaller than the spin healing length. We approximate the wavefunc- 
tion in the normal direction by a Gaussian with width d: ^m(^±, ^n) = '^m(^±)^(^n), where r± is the position vector 
in the 2D plane, is the coordinate in the normal direction, and h{rn) = exp[— /(4(i2)]/(27rd2)i/4. Multiplying the 
wavefunction to Eq. ([T]) and integrating over rn, we obtain the 2D GP equation. The equation is the same as Eq. ^ if 
one replaces ^777, with il^m^ as with r]as, Cdd with T^Cdd, and b with 6, where r] = J drnh^{rn)/ f drnh'^{rn) = l/V^Trd? 
and 

b, = I d\'^Q^^^\r^ - ri) [rM{F.)mnMr'A.)] , (24) 

with 

Qt^'^ir^ -r'^) = ]^IJ dr^dr'^h\r^)h\r'JQ^,{r - r'). (25) 
Starting from the 2D GP equation and following the above procedure, we derive the 2D hydrodynamic equation: 

+ (t^mass • V)/ = -/ X Beff, (26) 

with 

where i^mass and V are the two-dimensional vector and vector operator, respectively. When we consider a quasi-2D 
BEG, ntot, and i?mass are defined by means of iprn instead of ^rn- 



D. Dipole kernel 

This section provides the detailed form of the dipole kernel in 3D and quasi-2D systems under zero external field 
and under a strong magnetic field {p ^ Cdd^tot)- The derivations are given in Ref. |15j . 
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The dipole kernel in the laboratory frame of reference is given by 

Q(lab)(^)^ '^..^-3f,f. ^ (27) 

with r = |r| and r = r/r. The 2D dipole kernel in the laboratory frame is calculated by substituting Eq. (f27|) into 
Eq. ([25]) . For r±_ = (x, y) and rn = ^, the 2D dipole kernel is given by 

Q(f .'ab)(,^) ^J2e^,^--Q^2UM^\ (28) 

where 

4 / 1 \ f kl hky \ 

^(2D,lab) ^_ TT 10 +4^G(fcx^) , (29) 

^ \^ ^ -2) V 0-1/ 

with /c^ = (kx^ky), k± = kx,y = kx,y/k±^ and = 2/ce^^ e~^^ dt = y^/ce^^erfc(/c). It can be shown that 

G{k) is a monotonically increasing function that satisfies G(0) = and G{oo) = 1. 

When the linear Zeeman energy is much larger than the MDDI energy, we choose the rotating frame of reference 
in spin space by replacing with e~*^^*/^^^, and eliminate the linear Zeeman term from the GP equation. In 
this case, the contribution of the MDDI is time-averaged due to the Larmor precession, and we use the dipole kernel 
which is averaged over the Larmor precession period given by p3| 

1 1 _ 2 

Q^;^'Kr) = -- {5,^ - 35,,6,,) . (30) 

Substituting Eq. (|3Q|) into Eq. (j25j), we obtain the time-averaged 2D dipole kernel in the rotating frame as 

-^^^r^) = (V - 3(5,^,5,.) ^6^'=-'-- Qfe,, (31) 



where 



= y {l - 3(e„ • e^f - 3G{k^d) [(ef • k^f - (e„ • e'^f] } . (32) 
Here, k± = k±/k±, is the unit vector normal to the plane, and is the vector of projected onto the 2D plane. 



E. Stereographic projection 

The spin dynamics are now described by Eqs. ([20]) . ([2T]) . and ([23]) or ([26]) . We rewrite the equations by means of 
stereographic projection we employ a complex number (p = {fx + ify)/{l + fz) to express the spin variables, 

— . . 7 Jy — I 7 Jz — ^ • [o6) 



Equation ([20]) is rewritten as 



l + (/9(^* ' 1 + W* l + (/:^(^* 



V X i;^ass = ^F---^ ^, 34 

M (1 + (fif*)'^ 

while Eq. ([2T]) remains the same. In order to rewrite the equation of spins, we need to specify the dimensionality of 
the system and the direction and strength of the external field. In this paper, we consider the following two cases: (i) 
a quasi-2D system normal to the z axis under zero magnetic field, and (ii) a quasi-2D system normal to the y axis 
with a strong magnetic field along the z axis. Case (ii) corresponds to the situation in the Berkeley experiment [9[. 
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For case (i), we take = z and p = q = ^ and use Eqs. (|24|) , (|26|) , (|28|) , and (|29|) . Using the stereographic 
projection, the equation of motion of spins is given by 



\' ^ = -t^mass • + -— Vntot' V(^+—VV-ir7 1 ^ , 

at 2Mntot 2M M 1 + (^(^* 



2h 



^*(r') 



+ (/:)(r')(/:?*(r') 



2h 



h 

where the subscript ± was omitted for simphcity and 



o 

hi{k) = —-4TTG{kd), (36) 

h2{k) = 47rG{kd){k^^ikyf. (37) 

For case (ii), we take = = i, and p = 0, and use Eqs. (|24|) . ([26]) . (f3T|) . and ([32]) . Then, the equation of 

motion of spins is described as 

d(p(r,t) ^ ^ ^ ^^^2 ih (p*(\/(p)'^ 

^ = -^mass • + -— Vntof V(^ + —VV-T7 1 1 , 

at 2Mntot 2M M 1 + 



2ir]CddF 

-LP 



ig(2F - 1) 1 - w* 



^, (38) 



1 + ^' 

where ^ r = (x, z) and k± ^ k = {kx^ kz). 

III. DYNAMICAL INSTABILITY 

The hydrodynamic equations derived above give a rather straightforward approach to the analysis of the spin 
dynamics in a spinor BEC. In this section, we analyze the dynamical instability for cases (i) and (ii). Here we 
consider a uniform quasi-2D system and assume Vntot = 0. 

A. Case (i): Instability under zero external field 

Here we analyze the dynamical instability under zero external field for two initial stationary structures: uniform 
spin structures polarized normal to the x-y plane {(fo = 0) and in the x-y plane {cpo = 1). 

First, we consider the case in which the spins are polarized normal to the x-y plane, i.e., in the z direction, cpo = 0. 
Substituting cp = 0-\-S(p and Vmass = vq^Sv into Eq. (|35]) . we obtain linearized equations of 6(p and Scp"^. Performing 
Fourier expansions Sip = 5(pk^^^"^ and 5(p* = 5(p*_j^e^^"^ ^ we have 

d ( S(pk \ _ f -go - 9i -92 \ ( ^^k 



dt V ^r-u )~ n\ gl -90^91 )\ sr-j, J ' 
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-0.1 -0.05 0.05 0.1 

kx/27i (iim-'') 



FIG. 2: ReA+(/c) of the uniform spin structure polarized in the z direction under zero external field. The fluctuations in the 
black region are dynamically unstable and grow exponentially. Here, ntot is given by V^Trd^nsB with nsD — 2.3 x 10^"^ cm~^ 
and d = 1.0 /am. The other parameters are given by the typical values for a spin-1 ^^Rb atom: M = 1.44 x 10~^^ Kg, F = 1, 
gp = -1/2, and ^hf = 6.835 GHz x h. 



where 



9i{k) = ^ - 2^Cdd [2 - G(M)] , 



92{k) 



2M 

27r C(^(^G{kd){kx + iky)^ . 



Here, Cdd = ^Cdd^tot^ and k = {kx^ky). The eigenvalues of the 2x2 matrix in Eq. (|39|) are 



A±(/c) 



52 I 



9l 



(40) 
(41) 
(42) 

(43) 



The system becomes dynamically unstable when one of the eigenvalues has a positive real part; that is when 
ReA+(/c) > (or |^2p — > 0). The wavevector dependence of ReA+ is shown in Fig.[2l When the BEG is polarized 
perpendicular to the 2D plane, the MDDI is repulsive and isotropic in the 2D plane. Thus, the BEG is unstable 
against spin flip, and the unstable modes distribute isotropically in the momentum space. The unstable region in 
the momentum space has a ring shape. The radius and width of the ring are estimated as ko = {2/h)^/27TMcdd and 
Ak (v^/8)A:o<i(4 — ^/7^kod)^ respectively, for kd <^1. 

Next, we consider the uniform spin structure polarized in the x direction, cpo = 1. We obtain the linearized equations 
in a similar way to the above. Substituting (p = l-\-S(p and Vmass = Vq-\-6v into Eq. (|35]) and performing the Fourier 
expansion, we obtain the equation of the same form as Eq. (|39|) with 



9o{k) 
9i{k) 
92{k) 



hvo ■ k, 



2M 
2ncM 



27rcdd 



G{kd){l - kl) 



1 - G{kd){l + k^) 



(44) 
(45) 
(46) 



In this case, it can be shown that ~ 9i is always negative. Then, the eigenvalues which are given by the same form 
as Eq. (|43|) are purely imaginary regardless of k. Hence, the spin-polarized state along the 2D plane is stable under 
zero magnetic field. 



B. Case (ii): Instability under a strong magnetic field 



Here we analyze the dynamical instability for the helical spin structure, (fo 



^2(fc« r-u;«t)^ characterized by the 
helix wavevector Ka in the x-z plane under a strong magnetic field in the z direction. Substituting cp = cpQ and 
'^mass = '^0 into Eq. (|38|) gives uja = vq • K^a- Substituting (p = (po{l -\- S(p) and Vmass = vq -\- Sv into Eqs. ([2T 

and holding the terms up to the first order of the fluctuations, we have V • t;o = 0, V • = 0, V x t^o 
\/ X 6v = i^iySif + VSif*) X A^c^. After Fourier expansions of S^p., S^p* and 8v = Sv^e^^'^ , we have 



0, 



and 
and 



{k • K,(y)k 

P 



(47) 
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(a) uniform 



(d) helix 120)1171 




120 mG 



120 mG 



-0.04 -0.02 0.02 0.04 

kx/27i (iim-'') 



-(e) helix 120iim 



160 mG 



(b) uniform 



160 mG 



(f) helix 60iim 




120 mG 



(c) uniform 



200 mG 



(g) helix 60|im 







160 mG 



FIG. 3: ReA+(/c) for (a)-(c) uniform spin structures and (d)-(g) spin helices in the z direction with a pitch 2ti / under a 
strong magnetic field. The helical pitch for (d) and (e) is 120 /im, and that for (f) and (g) is 60 /im. The external field B is 
[(a), (d), and (f)] 120 mG, [(b), (e), and (g)] 160 mG, and (c) 200 mG. The other parameters are the same as those in Fig. (2] 



Substituting ip = ip^il^ Sep) into Eq. (|38|) and applying Eq. (|47|) , we obtain an equation of the same form as Eq. ([39|) 
with 



^o(^) 

gi{k) 

92{k) 



2M 2M 



{k • K'a)^ 



2M 



.2 

2 



2 

2 



Cdd 

4 



(2 



Cdd 

4 



qi2F - 1) 



q{2F-l) 



Qk- 



CddQk 



(48) 
(49) 

(50) 



Here, k = ikx^kz) and the Fourier transform of the dipole kernel is now simply given by Qk = (27r/3)[l — 
3{kz/kfG{kd)]. The eigenvalues of the 2x2 matrix in Eq. §9^ are given by Eq. (031) with Eqs. (|l8l)-(|S0l). When 
K,a//z, Vq = 0, and neither the MDDI nor the quadratic Zeeman effect exists (cdd = ^ = 0), the eigenvalues coincide 
with the dispersion relation derived in Ref. p^Tf . 

Figure [3] illustrates the wavenumber dependence of ReA+(A;) for uniform and helical spin structures under various 
magnetic field strengths, indicating the region of dynamically unstable modes. The dynamical instability discussed 
here agrees qualitatively with that obtained by the Bogoliubov analysis [HI, [25|. However, there is a quantitative 
discrepancy in the magnetic field dependence of unstable modes caused by the fact that the local magnetization of 
the condensate is assumed to be fully polarized in our method. However, when q is not sufficiently small compared 
with the ferromagnetic interaction, the amplitude of the magnetization decreases as q increases. For the parameters 
used in the calculation for Fig. [3l our assumption is valid for B <^ 480 mG. 



IV. MAGNETIC FLUCTUATION PREFERENCE 



We also investigate the magnetic fluctuation preference for two cases, that of dynamical instability under zero field 
for the uniform spin structure polarized normal to the x-y plane {cpo = 0), which is discussed in Sec. IIII A( and that 
of dynamical instability under a strong magnetic field for the helical spin structure {cpo = g<'*« '^-'^«^))^ which is 
discussed in Sec. IIIIBI 

For the case of zero external field for cpo = 0^ two kinds of magnetic fluctuations are considered: the x-direction 
fluctuation Sfx and the ^/-direction fluctuation Sfy. They are described by the first order of Sep: 5 fx = 2'Ke{5ip) and 
5fy = 21-m{5ip). Namely, the x- and ^/-direction fluctuations are characterized by the real and imaginary parts of 5ip^ 
respectively. 
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FIG. 4: (Color) (a) Magnetic fluctuation preference 0{k) for the unstable mode shown in Fig. [2] The x-direction fluctuation is 
dominant in the red regions (7r/4 < < 7r/2) and the ^-direction fluctuation is dominant in the blue regions (0 < < 7r/4). 
(b) (c) Schematic pictures of the magnetic patterns induced by the dynamical instability at points (b) P and (c) Q designated 
in (a). The arrows show the local magnetizations projected onto the x-y plane. The MDDI energy for configuration (c) is lower 
than that for (b). 



Let us reconsider the Fourier expansion of Scp^ 

= ^ [Ar sin{k ■ r + c^r) + iAi sin{k ■ r + ai)] + S(po, 



where 



(51) 



(52) 
(53) 



and ttR = tan ^[Iie{6tpk + S'fi-k)/^Txi{5(p-k — S(pk)] and ai = tan ^[lm.{5(pk + 5(p-k)/'Re{6(pk — S'P-k)]- As A+(fe 
0) = 0, Si^o = 0. We introduce the quantity 6, which characterizes the magnetic fluctuation preference: 



= tan-\AR/Ai). 



(54) 



Here, we calculate 9{k) for the unstable modes shown in Fig.|2j The eigenvector which corresponds to the eigenvalue 
A+ of the 2x2 matrix in Eq. is given by 



{9i+W\92\^-9i)/V^\ 



-91/VW 

where gi and g2 are defined by Eqs. (^T]) and (|^ . Then, from Eqs. (|52]) - (|55]) , we obtain 

1/2 

m 



(55) 



tan 



1^2 P - giRe{g2) + lm{g2)^/\g2\^ - gf 



I ^2 p + ^1 Re (^2 ) - Im {g2 ) \/|^2p -^i 



(56) 



We can consider the x-direction fluctuation to be dominant for 7r/4 < < 7r/2 and the ^/-direction fluctuation to be 
dominant for < 6> < 7r/4. 

The wavevector dependence of 0{k) for the unstable mode shown in Fig.[2]is illustrated in Fig.HJ^a). In the red (blue) 
regions, the x(^)-direction fluctuation is more dominant than the ^(x)-direction fluctuation. The schematic pictures 
of the magnetic patterns induced by the dynamical instability at points P and Q are illustrated in Figs. IH^b) and 
(c), respectively. At point P in Fig. |4l^a), x-direction fluctuation occurs and the wavevector of the magnetic pattern 
is directed in the x direction as shown in Fig. IH^b). The MDDI energy for this configuration is higher than that of 
Fig. m^c), where ^-direction fluctuation is induced. In other words, the MDDI energy is reduced by the magnetic 
fluctuation, and the reduction is larger in pattern (c) than in pattern (b). The reduction in MDDI energy is converted 
into kinetic energy (~ /c^), which is higher at point Q than at point P in Fig. HJ^a). This explains why the magnetic 
fluctuations change from x-direction to ^-direction along the kx axis. 
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200 mG 
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120 mG 



120 mG 



FIG. 5: (Color) Magnetic fluctuation preference 0{k) for uniform and helical spin structures in the z direction with a pitch 
I Ka /im under a strong magnetic field. The initial structure is uniform for (a) and (b), a helix with a pitch (c) 120 /xm, and 
(d) 60 /xm. The external field B is [(a), (c), and (d)] 120 mG and (b) 200 mG. The fluctuations are longitudinal in the red 
regions (7r/4 < Q < 7t/2) and transverse in the blue regions (0 < < 7r/4). The other parameters are the same as those in 
Fig. El 



Now, we investigate the magnetic fluctuations for the situation discussed in Sec. IIII Bl The fluctuations are con- 
sidered to be longitudinal or transverse. The longitudinal and transverse fluctuations are represented by Sfz and 
^{fx + ^/^y), respectively. Substituting (p = (po{l + Sep) and \(po\'^ = 1 into Eq. ([33|) . we obtain the expressions for the 
two types of fluctuations described by the first order of Sep: Sfz = —'Re{S(p) and S{fx -\-ify) = ^o^^{S(p). Namely, the 
longitudinal and transverse fluctuations are characterized by the real and imaginary parts of S(p^ respectively. 

We can also apply the above method to discuss the magnetic fluctuation preference, which is characterized by 
= tdin~^ (A-^/Ai)^ in the present case. Since g2 = we can simplify Eq. (|56]) as 

^(fe) =tan-\/ ^'~^\ (57) 

where gi and g2 are defined by Eqs. (|49|) and (|5Q]) . respectively. The magnetic fluctuation is longitudinal if 7r/4 < 
< 7t/2 and transverse if < ^ < 7r/4. 

The wavevector dependence of 0{k) for the unstable modes shown in Figs. [3] (a), (c), (d), and (f) are demonstrated 
in Fig. m When the dynamical instability has a round shape [Figs.[5fa), (c), and (d)], the fluctuations are transverse 
for small k and longitudinal for large k. Figure [5fb) looks more complex than the schematics for the other cases: 
fluctuations are transverse for k/ /z and longitudinal for k/ /x. The magnetic fluctuation preference is consistent with 
that obtained by the Bogoliubov analysis [l5|, [25| , and discrepancies appear in strong fields for the same reason as 
that for the dynamical instability. 



V. CONCLUSIONS 



Employing our hydrodynamic description derived in Sec. [Ill we have demonstrated some simple examples of the 
analysis of dynamical instability and magnetic fluctuation preference in Sec. IIIII and Sec. IIV| respectively. Once one 
finds a stationary solution of the hydrodynamic equations, it is a straightforward task to obtain the analytical form 
of the dynamical instability: the only necessary step is the diagonalization of a 2 x 2 matrix. The eigenvalues and 
eigenvectors of the matrix lead to the dynamical instability and the magnetic fluctuation preference, respectively. 
Although we have discussed just a few types of spin structures and external fields for simplicity, our method can be 
applied to other conformations. 

In conclusion, we have introduced the hydrodynamic equations for a ferromagnetic spinor dipolar BEG with an 
arbitrary spin by means of stereographic projection. This simple description provides a straightforward approach 
by which to investigate spin dynamics, i.e., dynamical instability and magnetization fluctuation preference, which 
are expressed in analytical forms. The description should also be useful for the study of the exact solutions of 
hydrodynamic equations of a spinor BEG. 
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Appendix A: Contributions from the short-range interaction, MDDI, and linear and quadratic Zeeman effects 



The contribution from the short-range interaction to the equation of motion of fz is calculated as 



dt 



1 

ih 



ih- 



2F 



dt 

s 



+ ih 



5—0, even 



Ms = -Slm'l' 



dt 

„ [-{ml\SMs){SMs\m'l')<i>i<i>*^,<i>^,^, 
+*:;(nl|5Ms)(5Ms|m'0*r*-'*H 



2F 



E E E"^("^^i^^«)(^^«i"^'^')[-*-*'*-'*^ + *-***™'*H 

5=0, even Ms^-Slm'l' 

Yl iy^s ~ m'){m, Ms - m\SMs){SMs\m, Ms - m)^l^^l^^_^^m^^Ms-m' 

Ms = — S mm' 



5=0, even 
2F 



E 



5=0, even 



0, 



27Th 

Jm 



as 



Z (m - mO(m, Ms - m\SMs){SMs\m' , Ms - m^K^k-m^-'^ 



Ms—m' 



Ms =—S ^mm' 



■ + l'){Ms - I, l\SMs){SMs\Ms - I', 1')^*m^_i^:^ Ms-i'^fi 



(Al) 



where we have used {Fz)mn = ^^mn- Here, [•••]s denotes that only those terms that come from the short-range 
interaction are extracted. In the following, this notation is applied to the contributions from the MDDI ([•••Jdd), 
linear and quadratic Zeeman effects. Since the short-range interaction (j2j) is invariant under spin 

rotation, [dfx/dt]s and [dfy/dt]^ also vanish, which are shown in a similar way by choosing the spin quantization axis 
along the x and y directions, respectively. 

The contribution from the MDDI to the equation of motion of spin is calculated as 



5/^ 



dt 



dd 



ih 



^ ^ rr. 
dt 



1 

ih 

Cdd 

ih 

Cdd, . f 
Cdd 



dt 



dd 



h 



(bx/)^. 



(A2) 



We have used the relations = F^ and [F^^F^y] = ie^j,\Fx. 

Suppose the magnetic field is applied parallel to the z axis. Then, the contribution from the linear Zeeman effect 
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is calculated as 



dt 



1 

ih 



mn 



ih 



^ ^ rr 
dt 



ih 



0%, 
dt 



P_. r 

P 



(A3) 



The contribution from the quadratic Zeeman effect is calculated in the same way as the that shown above. 



dt 



1 

ih 



mn 



ih 



^ ^ rr 

dt 



*1 ih 



dt 



ih 
ih 



^ Kf 



(A4) 



where M^y is a nematic tensor defined by 



(A5) 
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